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CN I. INTRODUCTION 

A new observational window to probe the structure of the universe is opening up: 21 cm redshift surveys. Such 
surveys offer the prospect of mapping out the distribution of neutral hydrogen in the universe, and thus will be able 
to probe the distribution of matter in the "dark ages" , before star formation sets in (see [JJ for an extensive review) . 
Regions in space which have an overdensity of baryons will lead to excess emission or absorption signals in 21 cm 
redshift maps. 

In a previous paper [5] we considered the signal of a single cosmic string in a 21 cm redshift map. Cosmic strings 
moving through space produce nonlinear matter overdensities in their wakes - even at redshifts larger than that of 
reionization. We pointed out that cosmic strings which are relevant to cosmological structure formation will leave 
a large imprint in 21 cm redshift surveys, in particular at redshifts larger than that corresponding to reionization, 
a redshift range where the "noise" from reionization processes is not present. As will be briefly reviewed below, we 
found that a single cosmic string gives a characteristic signal in 21 cm redshift maps: a wedge which is thin in redshift 
direction and extended in both angular directions (see Figure 1). The planar extent of the wedge is set by the length 
+3 of the string at the time it is formed, and the wedge thickness depends on the string tension. The angular extent of 
the signal from a string wake produced at the time of recombination is of the order of 1°, and the mean thickness 8z 
i— —i of the wedge in redshift direction (which depends on the tension of the string ) is about Sz/z ~ 10 -4 if the string 
tension /x is given by Gfi = HP 7 (G being Newton's gravitational constant) and the string is moving with a velocity 
close to the speed of light. The brightness temperature of the string signal depends both on the redshift of string 
formation and that of 21 cm emission, but only very mildly on the string tension. The 21 cm brightness temperature 
of a wake at a redshift z e + 1 = 20 can be of the order of —15 mK for strings formed during matter radiation equality 
with Gfi = 10~ 7 . Based on this large signal we expect that it will be much easier to see the signals of strings - if they 
C*") are present - in 21 cm maps than in other cosmological windows of observation. 

There has recently been renewed interest in the cosmological effects of cosmic strings (see e.g. [3H3 for reviews on 
cosmic strings and their effects in cosmology). It has been realized that many particle physics models beyond the 
Standard Model of particle physics give rise to cosmic strings. Such strings arise, for example, at the end of inflation 
t-H in many supergravity models 8 . Similarly, cosmic strings are a generic remnant of many models of brane inflation 
i- 

Furthermore, if strings exist, there will not be just one. In models which lead to cosmic strings, a network of strings 
will inevitably form in a phase transition in the early universe. Causality arguments |10) then imply that these strings 
will persist at all times up to the present time. Both analytical arguments [3 5j and numerical simulations [11 tell 
us that the network of cosmic strings will take on a "scaling solution" in which the average quantities describing the 
network are invariant in time if measured in Hubble length H (t), where H is the expansion rate of space. 

In this paper we will compute the angular power spectrum at fixed redshift produced by the distribution of cosmic 
strings emitting 21 cm radiation. Specifically, we are interested in the shape and overall amplitude of the signal. The 
paper is organized as follows: We begin in Section [H] with a brief review of cosmic string cosmology and observables 
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of cosmic string wakes in 21 cm experiments. In Section III we review the cosmic string toy model of 12 which we 



use, and construct the cosmic string wake profile and number statistics. In Section [TV] we calculate the 21 cm power 
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FIG. 1: Geometry of the 21 cm signal of a cosmic string wake. The left panel is a sketch of the geometry of the wake in 
space-time - vertical axis denoting time, the horizontal axis denoting one direction of space. The string segment producing the 
wake is born at time t, and travels in the direction of the arrow, ending at the position X2 at the time 2t%. The past light cone 
of the observer at the time to intersects the tip of the wake at the time S2, the back of the wake at time si. These times are in 
general different. Hence, the 21 cm radiation from different parts of the wake is observed at different red-shifts. The resulting 
angle-redshift signal of the string wake shown in the left panel is illustrated in the right panel, where the horizontal axis is the 
same spatial coordinate as in the left panel, but the vertical axis is the redshift of the 21 cm radiation signal. The wedge in 
21 cm has vanishing thickness at the tip of the wedge, and thickness given by 8v at the back side. 



spectrum for these cosmic string wakes. Section IV contains the main result of the paper. Readers who are familiar 
with cosmic string wakes or 21 cm observations can skip Section [TT] or Section |DT] and jump directly to Section |TV] In 
Section W\ we present our conclusions and put our work in context. 



II. COSMOLOGY OF COSMIC STRINGS AND THEIR WAKE'S BRIGHTNESS TEMPERATURE 

Cosmic strings are lines 1 in space with trapped energy. It is the gravitational effects of this energy which leads to the 
role which the strings play in the formation of cosmological structure, as first pointed out in [T3]. As relativistic objects 
they have tension equal to the mass per unit length and velocity v s of the order c. Their relativistic nature is also 
responsible for the characteristic string signatures in cosmological observations - the fact that space perpendicular to 
a long straight string is not flat, but conical with deficit angle a — 8ttG^l [Hj. This deficit angle leads to gravitational 
lensing of microwave radiation about the string which in turn leads to a line discontinuity in the temperature on the 
CMB [IS] . These line discontinuities can be searched for [16] using edge detection algorithms such as the Canny 
algorithm [T7] . Maps with good angular resolution such as those that are being produced by the South Pole Telescope 
(SPT) [TS] and Altacama Cosmology Telescope (ACT) [T^] are ideally suited to search for cosmic strings. 

A long straight cosmic string segment moving through space will lead to a wedge-shaped overdense region in its 
wake, a so-called "cosmic string wake" [20]. To understand what a wake is let us look at the physics from the point 
of view of an observer on the string. This observer will see matter streaming by on both sides. But since space is 
conical, the matter acquires a velocity 5v given by 

Sv = 4irG[iv s j(v s ) (1) 

towards the plane behind the string. This leads to the formation of a wedge of twice the background matter density 
behind the string. The wake, in turn, will grow in thickness by gravitational accretion, as studied e.g. in |21| using the 
Zel'dovich approximation. If a cosmic string is laid down at a redshift Zi, then the mean width of the wake (defined 
as the overdense region which has collapsed and undergone shock heating) at redshift z e will be 

w(t e ) = ^ G ^ s7 (^[j^^] 2 . (2) 
5 (z e + 1) 

The extra baryon density in the wake will lead to extra 21 cm absorption or emission. If the 21 cm radiation is 



These lines are in fact thin tubes with thickness given by a microphysical scale. 
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emitted from a wake at redshift z e its brightness temperature is given by 

5T b {z e ) = [70 mK]-^- (l - { \ + ■ (3) 

1 + x c \ Ik J 2 sin 8 

Here T 7 is the CMB temperature, Tk is the wake gas kinetic temperature, 8 is the angle of the 21 cm ray with 
respect to the vertical to the wake, and x c are the collision coefficients whose values can be obtained from the tables 
listed in [23] • Here and throughout we take the cosmological parameters to be H = 73 km s _1 Mpc -1 , = 0.0425, 
fi TO = 0.26. We work with a matter dominated universe for z < 3000 with the age of the universe to = 4.3 x 10 17 s. 

Note the Tk and the 8 dependence in Eq. [3] The string tension and the string's speed affect the temperature Tk 
of the neutral hydrogen in the wake, while the string's direction of motion determines 8. We simplify our calculations 
below by working with the brightness temperature of a string oriented so that 2 sin 2 8 = 1 and (v s j(v s )) 2 = 1/3. Thus 
the only string parameter that will determine the brightness temperature will be the string tension Gfi. 

Whether the string signal is an emission or absorption signal (positive or negative brightness temperature, respec- 
tively) depends on the ratio between the wake temperature of neutral hydrogen and the temperature of the microwave 
photons. The former is larger at lower redshifts, the latter is larger at higher redshifts. At redshifts below z = 20, 
reionization effects become important so we consider 21 cm radiation emitted at z e = 20 from a string formed at the 
redshift Zi + 1 = 3000 corresponding to the time t eq of equal matter and radiation. We found [2] that the critical 
value of Gfj, at which the transition from absorption to emission occurs is (G/x)e — 0.23, where (Gfi)^ is the string 
tension in units of 10~ 6 . 



III. COSMIC STRING WAKE NETWORKS 



A. Scaling solution 



The scaling solution is characterized by a random walk-like network of "infinite" strings (length larger than the 
Hubble radius) with step length comparable to the Hubble radius, plus a distribution of string loops with radii R 
smaller than the Hubble radius. According to numerical simulations [11) . the long strings dominate. For our study, 
we will make use of an analytical toy model for the distribution of strings first introduced in [12) : we divide the 
relevant time interval (in particular the time interval between t eq and the present time to) into Hubble time steps. At 
each step, space is divided up into boxes of Hubble size at the corresponding local times. In each time interval we lay 
down Nh straight cosmic string segments per Hubble volume. According to numerical simulations the value of Nh 
ranges from 1 to 10. Each string segment laid down at time t has a physical length at that time which is given by 
c\t (where c\ is a constant which is expected to be slightly smaller than 1, and whose value is expected to decrease 
as Nh increases) and depth given by v s j(v s )t, where v s is the velocity of the string in the plane perpendicular to its 
tangent vector, and "f(v s ) is the associated relativistic gamma factor. The value of v s is expected to be of order one 
in units where the speed of light c is set to one. 

The earliest Hubble time interval which we consider is at the time of equal matter and radiation, and the corre- 
sponding size of the spatial boxes at that time is ti = t(z oq ). The wakes produced by these strings are in fact both 
the most numerous (since a fixed comoving volume contains the largest number of wakes of that size) and the thickest 
(since the dark matter fluctuations have had the longest time to grow. The baryons are tightly coupled to radiation 
until the time i rcc of recombination but then fall into the potential wells created by the dark matter. String wakes 
laid down before t eq will have a thickness which scales as (z(ti) + l) -1 / 2 since the dark matter perturbations cannot 
grow while the dark matter forms the subdominant component of matter. 

Later time steps are centered at the values t^ = 2t(z eq ), . . ., t m = 2 m ~ 1 t(z cq ) respectively. Thus the total number 
of these Hubble time steps is 



M t , 



Int 



log 2 



(4) 



where Int(a;) is the integer part of x. Taking z cq ~ 3000, we have M to t = 17. However, as we discuss later in 
Section III E not all of these Hubble steps (at which the wakes are laid down) contribute to the 21 cm emission power 
spectrum. 

At each Hubble step, the cosmic strings are taken to be "created" at the beginning of the Hubble step and 
"annihilated" at the end of the step. In reality, the string network is a complicated dynamical distribution. Since 
the strings are moving relativistically, intercommutation events of two strings are frequent. For a fixed string, it lasts 
typically about one Hubble time step until this string intersects another one. Such an intersection will lead to a 
coarsening of the network of long strings via the production of cosmic string loops and the consequent straightening of 
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the remaining long string network. Hence, we arc modelling the distribution of the long string network in any Hubble 
time step as a collection of straight string segments which result from an initial intercommutation event at the time of 
"creation" , and undergo a next intercommutation at the time of "annihilation" . Since the intercommutation events 
are uncorrelated in space, it is a good approximation to consider the centers and directions of motion of the string 
segments to be randomly distributed, and uncorrelated from one time step to the next. Similarly, we expect a spread 
in the distribution of the velocities v s . 

The cosmic string wakes are produced by the motion of the strings through space. Once produced, the wakes will be 
approximately static in comoving coordinates, except for the growth in thickness due to gravitational accretion of the 
surrounding gas. We will consider the shape profile and number statistics of these wakes in the following subsections. 



B. Wake profile 

A cosmic string segment laid down at time ti will generate a wake whose physical dimensions at that time are 

c x t l x tiV s j s x AttGhUvsJs . (5) 

The dimensions ciU and tiV s ^f s span the two length dimensions of the wake and are independent of the string tension. 
They are both of the order of the instantaneous Hubble radius. The third dimension is the width of the wake, which 
is much smaller than the length because it is suppressed by the small parameter Gfj,. 

These three dimensions evolve with time. Firstly, the two length dimensions are comoving in space. Thus at a later 
time t e , their physical size is 

CiU (Jtt) x Vs1sU (fri) - (6) 

where 2$ = zfa) and z e — z(t e ). For simplicity, in some of the later calculations, we shall use l(zi,z e ) to denote 
the length of the wake at time t e if the time of formation is U, ignoring the order one difference between these two 
dimensions. 

The width of the wake will grow in time. According to linear perturbation theory if follows that the comoving 
width dimension grows linearly with the scale factor. We are interested in the width of the region which has collapsed 
and shock heated. The determination of this region can be done using the Zel'dovich approximation |22j and gives a 
result which is smaller than the naive linear perturbation theory result by approximately a factor of 4 [5] . The final 
result is given in pi). 



C. Wake Number Statistics: Method 1 



For later use, we will calculate the number density of cosmic string wakes on a fixed redshift hypersurface in two 
ways, making use of two different approximations. Both calculations agree to within a factor of order one. 

We shall use S 2 to denote the hypersurface with fixed redshift z e from which the 21 cm signals are emitted. The 
comoving radius r of this sphere can be calculated as 

r{t e )= r^-^ 3 ^ 3 -^ 3 ), (7) 

Jt c a ( t ) a 

where ao is the value of the scale factor at the present time to and can be set to 1 without loss of generality. 

At the time t m when wakes in the m'th Hubble time interval are laid down, the physical radius of the above S 2 is 

R m = a(t rn )r{t e ) = 3^ 3 (tJ /3 - # 3 ) . (8) 

Except at the final Hubble time steps close to the present time, this value of R m is much greater than the Hubble 
scale 1/H(t m ) = 3t m /2, which is also the typical length of straight string segments. Thus, for our counting procedure, 
the string segments can be treated as point particles. Note that the length of the string wake is of the same order 
as the Hubble scale. This breaks the particle approximation which is made here. However, as we show in the next 
section, using a full string word sheet calculation, the simple and intuitive approximation made here leads to the 
correct result, up to a constant of order one. 

The number of strings that pass through the S 2 per unit area per unit time is 



n s x v s x cos(9 s , 



(9) 
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FIG. 2: Illustration of how the number of relevant wakes is computed in Section [ill D| Here the z sw direction is chosen to be 
perpendicular to the string world sheet. And in the figure we have rotated the z direction to be the vertical direction only for 
the purpose of convenience. 



where n s , v s , 9 S denote the number density, averaged velocity and direction of motion of the cosmic strings, respec- 
tively. Thus, the number of strings that crosses the the S 2 is 

N S 2 = n s x v s x cos8 s x 47ri? 2 l x l/H(t m ) . (10) 

Since the cosmic string network achieves a scaling solution, it is more convenient to use number Nh of strings per 
Hubble volume instead of the number density of strings. Thus, we have 

iV S 2 = N H x v s x cos6» s x AttR 2 u x H{t m f . (11) 

This Ng-i is also the number of wakes on the sphere where the 21 cm signal with redshift z e were emitted. 
The two dimensional number density of cosmic string wakes at the emission time t e can be calculated as 

N S 2 _ 4N H v s cos8 s f z e + l \ 2 
1l2D - 4na(t e ) 2 r 2 ~ 9t 2 n [z m + lj ' ^ 

For our examples below we have averaged the direction speed of motion such that cos 2 S = 1/2 and v s = 1/2, and 
take Nh = 10, t m = 2.6 x 10 12 s for the first Hubble step (to = 1) corresponding to z = 3000. 



D. Wake Number Statistics: Method 2, a Worldsheet Calculation 

We calculate the probability that a string world sheet intersects the S 2 . The string world sheet is created at t m , 
and destroyed one Hubble time after that, with string length l 8 . Note that the length scale l s is of order t m . Thus 
the current method should be more precise than the previous one. 

For simplicity, we approximate the string world sheet as circular. The rectangular world sheet could also be 
calculable, but with the need of taking into account considerably more effects. We use r ws to denote the radius of the 
circular world sheet. To best mimic a rectangular world sheet, we take the area of the world sheet to be unchanged 
while deforming its shape. Thus we have 

^ws = c l*m X t m V s Js , r ws = tmy/ciVsls/w . (13) 

Consider one single circular cosmic string world sheet, as shown in Figure [2j Without loss of generality, we consider 
the world sheet inside a cylindrical volume. The cylinder has volume 2-kIj. Here we take L > R, where R is the 
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largest of the scales appearing in (13). The probability for the string world sheet (which is inside this clinder) to 
intersect the S 2 can be calculated as 

^3 J ^ dz ws dR ws ^ R ws d6 ws x 8(i? ws - RZ°' + r ws )8(iC° j + r ws - R ws ) . (14) 



where &(x) is the Heaviside step function, and R^°^ = \J R 2 m — z^ s . Keeping in mind that R m 3> r ws , the above 
integration yields 

(Probability) ~ . (15) 

L 6 

Thus, for N ws world sheets, the number of intersections is A^ ws 7r i?^r ws / L 3 . The number of world sheets 7V WS is related 
to the number density of world sheets (which equals to the number density of strings) by N ws = 27rL 3 n s . Thus, 
finally, we have 



N S 2 = 2ir 2 n s R 2 n r ws = (Ann^R^H- 1 ) ( 



y 3v s cos 9 S 



(16) 



When the factor y / 7rcT?v7^/(3t' s cos0 s ) is of order unity, this result is in agreement with the result (10). 



E. Cutting Off the Hubble Steps 



For the calculations in [2] to hold the wake must undergo shock heating. This occurs when the kinetic temperature 
inside the wake is at least 2.5 times greater than the background gas temperature (the temperature of the matter 
particles outside of the wake, not the photon temperature T 7 ). Shock heating is more likely to have occurred for 
strings laid down earlier. Below we determine the upper cutoff time for shock heating to have occurred. 

Note that the later the wake is laid down, the lower the kinetic temperature. This leads to a lower bound on the 
redshift z m , after which the cosmic string wakes will not have undergone shock heating. The kinetic temperature for 
the gas inside the wake is [2] 



Until rcionization, the background gas is cooling adiabatically below redshifts of z s» 150 as: 



T g = 0.02K(l + z e ) 2 



The requirement Tk > 2.5T„ leads to 



1 > 



(1 + ^e) 3 



m(Gfi) 2 (v sls y 



(17) 



(18) 



(19) 



Only the Hubble steps with z m satisfying Eq. ( 19 1 lead to wakes where shock heating occurs 



This requirement that only shock heated wakes contribute significantly leads to a maximal number of Hubble steps 
in the toy model. The Hubble steps that should be taken into account are m = 1, . . . , M, where 



M = Int 



log 2 



= Int < log 2 



+ 1 



3/2 



(20) 



where Int(x) denotes the integer part of x. Later wakes which have not undergone shock heating are thinner (see 
Eq. [2]) and less numerous (see Eq. 12) than the earlier ones. In addition, they are not gravitationally bound due to 
free streaming. For both of these reasons, their contribution to our power spectrum is negligible. 

For example, when we take the 21 cm emission redshift to be z e = 20, and set the maximal redshift of cosmic string 
wake laid down time to be the matter radiation equality time z™ ax ~ 3000, then the number of Hubble steps which 
lead to important cosmic string wake 21 cm emission is as plotted in Fig. [3j 
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FIG. 3: The number of Hubble steps, for a given string tension, during which a wake emits 21 cm radiation. The parameters 
used are z e = 20, z™ ax = 3000 and v'i'yl = 1/3. 

IV. POWER SPECTRUM 

We will compute the angular power spectrum of the 21 cm brightness temperature at a fixed redshift with a string 
tension larger than the critical value for emission. A few words of warning are in order. Firstly, for a Gaussian 
distribution the two point function contains all the information about the distribution. However, in the case of cosmic 
strings, the signal is very non-Gaussian. Hence, the two point correlation function will miss important information. 
Secondly, one of the advantages of 21 cm surveys compared to other measures (e.g. CMB temperature maps) is that 
the maps are three dimensional and thus can contain much more information. In restricting attention to a fixed 
redshift, we are not exploiting this advantage. Nevertheless, in order to compare the 21 cm signal of cosmic strings 
with that of other sources, it is useful to compute the angular correlation function. We plan to extend our analysis 
and compute the three dimensional correlation function in upcoming work. 

We are performing our calculation in the framework of the toy model distribution of long strings of [12] reviewed 
earlier. In this analysis, the strings at different Hubble time intervals are statistically independent. Hence, we can 
first consider a fixed string formation time Hubble interval and sum the results over all redshift intervals which make a 
contribution. In the following subsection, we therefore consider the contribution from strings in a fixed single Hubble 
string formation time step. 

A. A single Hubble step 

First we calculate the two point correlation function (5Tb(0)5Ti 1 (x.)} The brackets indicate an ensemble average. 
We assume the validity of the ergodic hypothesis so that ensemble averages can be computed by volume averages. 
We shall use the flat sky approximation because we are interested in scales of order ten degrees on the sky or smaller. 
The position x in the two dimensional sky is the 2-d angular coordinate times the radial distance of the fixed redshift 
surface z e under consideration. 

Whereas in the three dimensional redshift-angle space a wake is long in two dimensions and narrow in one, the 
intersection of the wake with the fixed redshift surface produces an object which is long in only one dimension and 
typically narrow in the second. We will call the size of this second dimension the "projected width" . As illustrated 
in Fig. [4j the projected width w of the wake can be written as 

_ _ Jw/cos-0 for < i> < ir/2 — arcsin(w//) , . 

W ~ \ /for tt/2- arcsin(w/0 < V < tt/ 2 ' ^ ' 

where ip is the angle of the wake relative to the fixed redshift surface. 
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Since the positions and orientations of different string wakes at a fixed Hubble time step are taken to be random, 
the computation reduces to the computation of the ensemble average of the contribution of a single string wake. We 
will consider the ensemble average of the contribution of a single string in a box of dimensions L x L , where L is lager 
than the length of the wake we are interested in. Since the dominant contribution to the 21 cm signal comes from 
strings created at about i cq , we do not need L to correspond to a large angular scale. As in the studies of position 
space signals of cosmic strings to CMB temperature and polarization maps, an angular extent of 10° is more than 
enough (larger areas will obviously lead to reduced statistical error bars). Returning to the study of the effect of a 
single string wake, the geometry is illustrated in Fi g. |4| H ere we assume that the cosmic string wake intersects the 
surface of fixed redshift z e . As discussed in Section |III B| the number density of such string wakes is ti-id- We now 
study the contribution of this string to the correlation function. 

When both points and x are on the wake, the product JT^O^T^x) contributes ST b to the ensemble average. 
Otherwise, when either points or x are not on the wake, the product is zero. In this way, the calculation of the 
ensemble average is reduced to the calculation of probability that both points and x are on the wake. We have 

(6T b (0)6T b (x)) = Iff x (Probability) . (22) 

There are probabilities for both the horizontal and vertical coordinates of the points and x to be on the wake. When 
< 4> < arctan(Z/iu) (0 is labeled in Fig. |4|, the probability we calculate is a product of these two probabilities: 

I- w/ cos <f,-r I cos 4> f ~ I , 

(Probability) = Z X ^ *° r r < W J C0S< P (23) 

v I for r > uj/cos</> v ' 

where r is the magnitude of x and w = w/ cos ip is the projected width of the wake onto the fixed redshift hypersurface. 
The angles ip and </> are shown in Fig. [4] They should be integrated over when performing the ensemble average. 

On the other hand, there is the other region axct&n(l / w) < <f> < 7r/2, which is small on small scales but will give a 
comparable contribution on large scales (scales of order the length of the wake). In this region, the probability that 
both points and x are on the wake can be written as 

f j sin <p-r w/ sin f , ■ i 

(Probability) = x ~ r~ [ 0T r < sm J . (24) 

v for r > I sm o 




r 



LxL box on the fixed z e hypersurface 



FIG. 4: This figure illustrates the cosmic string wake profile. Both panels represent cross sections through space at the time 
t e . The "fixed redshift hypersurface" is the intersection of our past light cone with this space. The left panel introduces the 
angles which describe the relative geometry between the wake and the fixed redshift hypersurface, on which the 21 cm signal 
is emitted. The right panel shows the projection of the wake onto the fixed redshift surface. Here the "width" in the figure is 
w, the projected width of the wake. 



9 



Based on rotational symmetry of the ensemble, the correlation function will only depend on r, the magnitude of x. 
For a single wake, the considerations of (22), (23) and (24) yield 



(ST b (0)ST b (r)) 



26T b 
ttL 2 

26T1 2 
irL 2 



tt/2 



dip 



r/2 



dil> 



arctan(Z /w) 

tt/2 
arctan(Z / w) 



<fi I sin ip (w — r cos <fi)0(w/ cos 4> — r) 
ws,m.ip{l — r/sin</>)0(Zsin0 — r) . 



(25) 



where 0(x) is the Heaviside step function. Note that we have put a normalization factor 2/ir in front of the integration, 
because we are averaging over 1/8 of the solid angle by restricting ip and <p to both run from to tt/2. 

The angular power spectrum of zero mean fluctuations is obtained by taking the Fourier transform of the above 
correlation function and subtracting the square of the brightness temperature times (27r) 2 <$|>(fc) for the zero mode. 
Hence for k ^ the contribution of a single string wake can be written as 



/>oo t*2ir 

Pi(fc) = / dr dO 
Ja Jo 



ikr c 



} (ST b (0)ST b (r)) 



where r, 9 are the polar coordinates of x. One can insert Eq. ( 25 ) into Eq. ( 26 1 , integrate out 



(26) 



and ip to obtain 



the power spectrum on length scales much smaller than the wake length 
be well approximated by 



In this regime, the power spectrum can 



Al5T b 
L 2 k 3 



[1 — cos(wfc) — wkSi(wk) + wkSi(lk)] , 



(27) 



where Si(z) is the sin integral defined as Si(z) = Jq dt (smt)/t. On the other hand, on scales larger than the wake 
length I, the power spectrum P\(k) is independent on k. This can be seen by taking the limit k — ¥ in Equation (26 1 . 
The full power spectrum can be calculated numerically. 

Now we will consider multiple cosmic string wakes which are laid down within the same Hubble time step. Up to a 
factor of order one, these cosmic string wakes have the same dimensions as calculated in Section |III B| If we assume 
these wakes do not have any cross correlations between each other, then we simply have to multiply the result for 
a single string wake by the number of string wakes in the observation area, which is L 2 n2D- Thus, the result for 
multiple string wakes (single Hubble time step) in the region kl » 1 becomes 



P{k) 



4ln 2D ST b 
k 3 



[1 — cos(u>fc) — wkSi(wk) + wkSi(lk)] 



On scales much smaller than the wake length I, the dimensionless power spectrum can be written as 

-2 



A(k) 



2^ 



P(k) 



2ln 2 D&T b 
nk 



[1 — cos{wk) — wkSi(wk) + wkSi(lk)] 



(28) 



(29) 



Eq. (29) and the corresponding numerical result (shown in Fig. [5]) are the main result of our paper. 

Since, as mentioned earlier, P(k) is constant in the limit of small values of k, A(fc) scales k 2 for small k. The scaling 
for large values of k is fc" 1 , which can be obtained by expanding Eq. (29) in a power series. 

It is also convenient to use wk as the argument and rewrite the dimensionless power spectrum (for kl > 1) as 



A(fc) 



i-c»(-)- M *si(»* )+ -si|A*„* 



(30) 



Eq. (30) is plotted in Fig. [5] Note that the dimensionless power spectrum is almost scale invariant between the 
scales 1// < k < 1/w corresponding to the wake length and the wake width. When k < 1/i it decays as k 2 , and when 
k > 1/w it decreases as fc" 1 . 



2 Here dimensionless means that there is no k dimension. The power spectrum still has dimension of temperature squared. 
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FIG. 5: The dimensionless power spectrum A(fe). We plot A(fc) / \wln2D&Tb ) as a function of wk. In the plot we take l/w to 
be 9000. This value is calculated using z e = 20, z m = 3000, (G/x) 6 = 0.2, v s = 1/2, and cos 2 6 S = 1/2. 

As an example, we consider (G[i)q = 0.2 (around the current upper bound on the string tension), z e = 20 and 
Nh = 10 strings per Hubble volume. From Fig. [3] we see that only wakes from a single Hubble time step undergo 
shock heating. The relevant length scales of the wakes are 



w = 2.3 x 10~ 4 Mpc , I = 2.1 Mpc . 
The number density of wakes from Eq. |12[ and the induced 21 cm brightness temperature from Eq. [3] are 



(31) 



(32) 



n 2D = 0.12 Mpc~ 2 , ST b = -15 mK . 

In this case, the dimensionless power spectrum becomes 

A(fc) ~ 8000^K 2 for l/l < k < l/w . (33) 

Comparing this result with the square of the local brightness temperature from ^ we see that the power spectrum 
is suppressed by a significant factor. This is a reflection of the highly non-Gaussian distribution of wakes in space, 
and the resulting highly non-Gaussian distribution of the string-induced wedges in the 21cm brightness maps. The 
fraction of the 21cm sky at a fixed redshift which is covered by string-induced wedges is proportional to w/l ~ G/i. 
Thus, whereas the local brightness temperature inside a 21cm wedge is to a first approximation independent of the 
cosmic string tension, the two-dimensional power spectrum decreases in amplitude as Gfi decreases. The dependence 
on G/i is linear as can be seen from ( 30 ) , since in the flat portion of the spectrum the last term in the square brackets 
dominates. 

Finally we also would like to comment which features of the power spectrum are expected to exist for a realistic 
cosmic string network and which features are artifacts of the toy model. 

The decay rates of A (A;) at small and large values of k should be robust because they do not depend on special 
assumptions made in the cosmic string toy model. We also trust the nearly scale invariant behavior between 1 jl and 
l/w. 

However, the oscillations seen in the spectrum near k — l/l are probably an artifact of the toy model. This is 
because we have assumed all string wakes are laid down at one particular redshift z m . This is not accurate (up to 
an order one factor) even in the case that only one Hubble step is needed. String wakes laid down at different z m 
are expected to lead to cancellations because of the different phases of oscillation. We also expect that the spectrum 
should have a slight blue tilt near k = l/l. This is because the larger scale wakes are laid down at later times and 
thus lead to a lower brightness temperature. 



B. Multiple Hubble steps 



We now consider the case when there are several Hubble steps which lead to wakes with Tk < 2.5T g . In this 
situation, the power spectrum becomes a superposition of the power spectra in different Hubble steps. For example, 



11 



when we are considering measurements for emission redshift z e — 20, then when (Gfi)e > 0.22, multiple Hubble steps 
are necessary. 

We assume no cross correlation between string wakes in different Hubble steps. Then, the contribution from different 
Hubble steps can be added together directly. In the left and right panels of Fig. [6j we plot the dimensionless power 
spectrum from wakes at each Hubble step and the sum, respectively. Note that the first Hubble step has the largest 
brightness temperature. Thus the first step dominates the total power spectrum. 




k/Mpc" 1 k/Mpc" 1 

FIG. 6: The dimensionless power spectrum from string wakes at individual Hubble steps (left), and the sum over the steps 
(right). The parameters are z e = 20, 2™ ax = 3000, (G/x) 6 = 0.3, v s = 1/2, cos 2 6 S = 1/2. and N H = 10. 



V. CONCLUSION 

We have computed the angular power spectrum of 21 cm emission (at a fixed redshift) from a scaling distribution 
of string wakes. Wakes from strings present near the time of equal matter and radiation, dominate the spectrum. 
The dimensionless power spectrum scales as k 2 on long distance scales and as on short distances. In the above, 
"long" is relative to the length of the string wake which corresponds to an angular scale of about 0.1° and which is 
independent of the string tension. "Short" is relative to the width of the wake which depends on the string tension (its 
order of magnitude is G/i times the length of the wake. On intermediate scales, the spectrum is roughly scale-invariant. 

The above features of the 21 cm power spectrum from string wakes are robust against details of the string wake 
scaling solution. The overall amplitude, however, depends quite sensitively on the details. As an example, let us 
consider the redshift z e = 20, larger than the redshift of reionization so that we do not need to worry about "noise" 
from reionization processes. Then, for {G^)q = 0.2 (the current upper bound on the string tension, and assuming 
10 string segments per Hubble volume, the amplitude of the dimensionless power spectrum over the flat part of 
the spectrum is of the order A(fc) 1 / 2 ~ 90 /iK. This amplitude is significantly smaller than the local brightness 
temperature in a string-induced 21cm wedge in position space. Whereas the local brightness temperature in position 
space is to a first approximation independent of the string tension, the two-dimensional power spectrum decreases 
linearly with Gfi as the string tension decreases. This difference between the local and the root mean square power 
of the signal is a reflection of the highly non-Gaussian distribution of cosmic string wakes, and it reinforces the lesson 
learned from studies of cosmic string imprints on CMB temperature maps - namely that it is easier to detect the 
signatures of cosmic strings when analyzing the data in position space. 

It will be of great interest to study in detail the detectability of this signal with upcoming and planned experiments. 
To put the predicted amplitude into some context, let us recall that the current strongest upper bound on a cosmolog- 
ical 21 cm signal is 70 mK [S3] on angular scales of about 3.3°. The redshift range probed in this experiment (GMRT) 
lies between 8.1 and 9.2. One foreground for our signal is the power spectrum induced by mini-halos, nonlinearities 
which result from the Gaussian primordial fluctuations. An analytical estimate of the effect was given in , with 



12 



the predicted spectrum A(fc) 1 / 2 at redshift z — 15 peaking at k ~ 10 3 /iMpc^ 1 at an amplitude of about 100 mK and 
decaying linearly with k for smaller values of k. A detailed numerical study of the expected signal of reionization 
effects at a redshift of z = 11 was performed in |27j yielding a spectrum which peaks at lOmK at a value of the angular 
quantum number I of I ~ 10 4 and decaying roughly linearly with I for smaller values of I, and turning flat for larger 
values. Turning to low reshift surveys, we recall that 21 cm experiments are planned to be able to pick out a 20/iK 
signal from baryon acoustic oscillations from a noise of 300mK from large-scale stucture on a length scale of about 
8Mpc [28 . Thus, it appears to us that the predicted spectrum from a scaling distribution of cosmic string wakes is 
large and has a good chance of being visible in the planned high redshift surveys. 

Our work is based on a simple toy model for the distribution of cosmic strings which captures the main properties 
of string networks. Cosmic string loops and wiggles on long strings are, however, not contained in the present toy 
model. It would be of interest to extend our work to capture these aspects of real string networks. 

The restriction to fixed emission redshift cuts down a lot of the information which is contained in 21 cm redshift 
maps. Hence, it will be of great interest to extend our analysis and calculate the three dimensional redshift power 
spectrum of 21 cm emission from cosmic strings. Work on this topic is ongoing. 
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